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1. Introduction 


The new observations of the heavy ion collider facilities represent a great challenge to particle 
field theory. The achieved volume and life time of the produced hot and dense matter may now 
enable thermalization. On the other hand, due to the developments on the theory side, calculations 
of nonequilibrium processes become increasingly feasible. 

It has been recently observed that the phenomenological description of particle flow based 
on ideal hydrodynamics is extremely successful at RHIC energies |[T]]. The used models assume 
local equilibrium already before 1 fm/c. The competence of ideal hydrodynamics indicates a very 
early thermalization of the produced plasma 0. A further signal of an at least partial equilibration 
is that the global abundances from collision experiments obey simple statistical models assuming 
chemical equilibrium [^. 

The observed early onset of equilibrium calls for a novel theoretical paradigm. The perturba¬ 
tive thermalization time estimates are far beyond the experimental expectations. As a first step one 
has to understand, what level of equilibration is required to interpret the experimental data. The 
success of ideal hydrodynamical description might be explained by the early onset of an isotropic 
equation of state |Q] that may be present even in a far-from-equilibrium quantum field iJH]. Then one 
seeks for explosive processes that drive the system towards isotropy and equilibrium. A promising 
candidate is the development of plasma instabilities that may dramatically accelerate the evolution 
of anisotropic fields 

Nonequilibrium methods are also asked for in early Universe scenarios. Much of the efforts 
have been put into the quantitative understanding of the theory of reheating, which explains the 
rapid re-population of the dilute Universe after inflation [^. The final temperature is an important 
input to models of baryogenesis |^]. 

The intriguing questions of non-thermal field evolution cannot be answered by a simple per¬ 
turbative analysis. Even in a weakly coupled theory high orders become relevant with the evolution 
of the time. The elapsed time appears next to the coupling constant, making the effective coupling 
arbitrarily big at late times. This phenomenon, called secularity, invalidates approximations that 
truncate higher order diagrams in the equation of motion of the Green’s function |[T^. A way to 
evade this trap is to use self-consistent schemes 0]- 

There has been important progress in our understanding of nonequilibrium quantum fields 
using suitable resummation techniques based on two-particles-irreducible (2PI) generating func¬ 


tionals [12]. It has been shown by Cox and Berges [13] that the Kadanoff-Baym equations are 


suitable for direct numerical treatment. They used the 2PI effective action formalism [14] to calcu¬ 
late the self energies in as systematic way. They found that the studied scalar field reaches thermal 
equilibrium even in l-i-l dimensions, where binary collision are forbidden by kinematics on the 
level of the Boltzmann equation. 

The 2PI effective action techniques (or <I> derivable approximations) are known to provide 
time reversal symmetric, energy conserving equations for the field propagator ][TT[]. The equations 
involve the resummation of those diagrams that carry the time in the effective coupling, this cures 
the secularity of the perturbative treatment. 

The 2PI equations of motion have been solved for various models since the pioneering work 
by Cox and Berges [13]: ranging from the simplest <I>^ scalar theory [32, 29, ^ through the scalar 
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Figure 1: Evolution of the particle distribution function in the chiral quark model pl|]. To the left, the 
fermionic particle spectrum is shown for various times. The spectra have been transformed so that a linear 
corresponds to the Fermi-Dirac statistics. To the right, the final distribution is displayed without the transfor¬ 
mation. The coupled system of fermions and scalars dynamically realizes the Fermi-Dirac and Bose-Einstein 
statistics. 


0{N) theory with next-to-leading order large N resummation [I5,I^n,18,]^,^to the chiral 
meson model [21], often considered as a prototype of QCD. 

The 2PI formalism was used to show for the first time from first principles the formation of 
the Bose-Einstein and Fermi-Dirac statistics in a system of coupled fermions and scalars [^]. 
The used equations do not involve any of the statistical factors or constraints for the occupation 
numbers, not even the concept of particles. The resulting distributions are purely an outcome of 
the coupled field dynamics (See Fig. [I]). 

In fhe nexf sections we firsl review fhe 2PI effective acfion formalism on fhe example of a 
scalar fheory. Then we discuss fhe various resummafions necessary fo calculafe fhe self-consisfenf 
n-poinf functions fo fhe desired accuracy. The 2PI resulfs have been checked bofh in equilibrium 
and ouf of equilibrium, we briefly menfion fhe main resulfs of fhese fesfs. We also highlighf some 
of fhe numerical resulfs on fhe basic time scales of nonequilibrium evolution. 


2. The formalism 


The objecfive is fo give an equafion for fhe propagafor and also for fhe higher n-poinf funcfions 
of an evolving field. The 2PI formalism does nol use fhe concepf of quasi-particles, if esfablishes 
fhe relafionship befween various n-poinf funcfions in a sysfemafic fruncafion scheme. To keep fhis 
review simple, we resfricf our discussion fo fhe scalar propagafor, which already gives accounf for 
fhe evolufion of parficle disfribufion. The equations are worked ouf for higher n-poinf funcfions 
in Ref. [22]. The 2PI formalism for multiple scalar fields is discussed in defail in Refs. [16, S], 


for fermionic fields if has been elaborafed in Refs. Q. For renormalizable fheories fhe 2PI 
equafions of mofion are also renormalizable, fhe correcf way fo define a finife fheory has been 


discussed by many aufhors [25, 26, 22, 23, 24]. The inferesfed reader is referred fo fhe in-depfh 


review in Ref. [12]. 
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The proper two-point function, that we are looking for, is defined in terms of the one-particle 
irreducible (IPI) effective action: 


^-\x,y) 


5^{x)5^{y) 


( 2 . 1 ) 


In the case of more than one field components is a matrix in the field indices as well as in the 
space-time coordinates. 

The r[0] effective action is now a functional of the space and time dependent background field. 
A straightforward approximation would be a coupling expansion in r[0]. A danger in this approach 
is, that a Taylor expansion in (j) simply cancels the higher proper n-point functions. Negative results 
warn us that it is required to include high powers of (j) in the effective action. This can be achieved 
by imposing the truncation to the Legendre transform of r[0]. This way we selectively resum IPI 
diagrams to infinite order. 

The two-particle irreducible (2PI) effective action is defined as a Legendre transform of the ef¬ 
fective action with respect to the source .S'(x,y) of the two-point composite field operator {^{x)^{y)): 


r2Pi[<l),n<l>,K]] =r[0]-iy‘d4x Jd^ym(t>,K]i^,y)K{x,y) + (t>{x)K{x,y)(t>{y)) . (2.2) 


The subscript refers to the closed time path (CTP) contour for the time integral. 

The second argument of the 2PI effective action r2pi[0,G] is a generic two-point function 
G(x,y). For any given (j) background and vanishing sources 


^r2Pi [ip,G] 
5G{x,y) 


(2.3) 


with G = '^[(p]. From r2pi it is straightforward to go back to the IPI effective action using the 
identity r[0] = r2Pi[0,^[0]]. 

Unfortunately, the exact r2Pi functional is not known. Following Cornwall, Jackiw and Tombu- 
lis |[U one can decompose the 2PI effective action into tree-level (1st term), one-loop level (2nd 
and 3rd term) and higher order (4th term) contributions. 


r2Pi [0,G] = S[(/>]-h-tr^ [log [G ^]]-h-tr^ [Gq ^[(/»]G] (2.4) 

-|-r2 [<^,G] + const . 


Here Go[(/>] denotes the free propagator on the background, and S is the classical action. The decom¬ 
position (^^ defines a term r2[^,G] that can be entirely represented by two-particle irreducible 
diagrams, i.e. these diagrams do not fall apart if any two of its lines are cut. 

We will use a slightly modified decomposition, where the background dependence of the free 
propagator is put into the interaction part rint[0,G]: 

r2Pi [(/),G] = S[(^] + ^tr^ [log [G-i]] + ^tr^ [Gq i[(/) = 0]G] (2.5) 

d-Fint [^,G]+const . 

In the 2PI formalism the truncation is realized by the formal expansion of rint[^,G] in the 
coupling constant or in an other small parameter, such as the inverse number of field components. 
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There are Feynman rules to construct r,„f from the selected diagrams. The most important ingre¬ 
dient of these rules is to write a G{x,y) function for each propagator. For the vertices, however, the 
bare couplings appear (with counterterms, eventually). 

For the most trivial case of scalar A<I>^/24 theory one simply keeps: 

rint[G] = —J dx^{x)G{x,x)^{x) — — J dxG^{x,x) 

Jdxdy(^{x)G^{x,y)(j){y) + i^ J JdxdyG‘^{x,y). (2.6) 


The first two terms correspond to the (collisionless) Hartree approximation [^], the next two terms 
give account for the “particle collisions” at lowest order. An expansion to higher orders in A is 
found in Ref. |^]. 

When using a truncated r 2 pi we have to rethink the exact identities derived so far. We define 
fhe resummed two-point function G[(/»](v,y) as the solution of Eq. ( ^ ) with the truncated r 2 pi. 
Using this as an argument of r 2 pi one uses 


r[(/)]=r2Pi[(/),G[(/)]] (2.7) 

as a definition for the truncated IPI effective action pT|]. As already noted, this relation turns to an 
identity in the case without truncation. This truncated r[0] is used then to define W according to 
Eq. (^TTI). This ^{x,y) propagator is in general, not identical to G[(/)](x,y) 

As a first step we have to solve Eq. (^^, which can be put to the form of a Schwinger-Dyson 
equation using Eq. ( |2.5[ ): 

G \v,y) =Go^[0 =0](v,y)-r[G[^]] (v,y) (2.8) 


with 


(-^a) = 2/ 


stm m 

5G{x,y) 


(2.9) 


Here the self energy (Z) is built of one-particle irreducible diagrams constructed form resummed 
lines (G[0]) and bare vertices. (Coming from the first term in Eq. ( ^^ it may actually contain a 
diagram ~ 5^{x,y)(j){x)(j){y), which is not IPI. Using the original decomposition (^^ this term 
would appear in G^^ [0].) 

In the second step we evaluate the 2PI effective action at G[(p] and calculate as follows: 


^ ‘[^](v,y) 


■ 

5^{x)5(j){y) 


. ^^r2Pi[0,G[(/)]] 
' 5(j){x)5(^{y) 


( 2 . 10 ) 


Analogously, any higher n-point functions are available if G[(j)] is known for any (j) background. 

In the symmetric phase (0 = 0) Eq. ( 2.10| ) takes a simpler form if the model has the Z 2 {(p < 
—(p) symmetry: 


'[0 =0](v,y) =Gq^[(P =0]{x,y)-i 


52ri„t[0,G[(/)=O]] 


5p{x)5p{y) 


( 2 . 11 ) 
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We emphasize at this point, that in the broken phase additional terms appear 
additional terms the propagator 


5 . Without these 


would violate Goldstone theorem [ 131| , |16| ]. 

There are some “friendly” truncations (including Eq. (|^) for which Eim has the symmetry: 


5^rint[0,G] 


= 2 


5^rM[<P,G] 


l )=0 


5G{x,y) 


d^{x)5^{y) 

This further simplifies Eq. ( 2.11| ) to 

^[(/) = 0](;c,y) = G[(^ = 0]{x,y). 


i=0 


( 2 . 12 ) 


(2.13) 


Eor a Z 2 symmetric exact theory Eq. (|2.12| ) is always granted, as expected, and G are then 
equivalent. 

All the models numerically investigated so far (the three-loop <I>^ theory, the 0(A) model 
truncated at next-to-leading order in 1/A and the two-loop chiral quark model) have this nice 
feature. This explains why numerical works directly use the solution of the Schwinger-Dyson 
equation 1 


3. Evolution equations 


Eq. ( |2.8[ ) uses complex and nonanalytic functions on CTP contour, thus, it is not suitable for 
a numerical treatment. Therefore we introduce the statistical propagator F{x,y) and the spectral 
function p{x,y) and the corresponding functions for the self energy, too [32]: 


G{x,y) = F{x,y)-^p{x,y)sgn^{xo,yo), (3.1) 

^[G](^,y) = -i^o{x)d<^{x,y)+I.^{x,y)-^LP{x,y)sgn^{xo,yo)- (3.2) 

The real (F) and imaginary (p) part capture different physical information: p gives what states are 
available and how stable they are, while F tells how they are occupied. In equilibrium, the F and p 
functions are connected by the KMS condition, but they are a priori independent out of equilibrium. 
These functions do not depend on contour variables, all discontinuities are encapsulated in the 
contour delta 5^ and sign function sgn^. One can transform Eq. ( ^^ into integral equations, 
where the integral is defined on fhe closed-fime-pafh contour [^]. Using fhe F and p functions in 
fhese infegral equations fhe disconfinuifies appear as boundaries of fhe fime infegral. Eor a scalar 
field fheory, fhe equations of motion in terms of fhe F and p funcfions read: 


{d^ + m^ + LQ{x))F{x,y) = j dzL^{x,z)p{z,y) - j dzlP {x,z)F{z,y), (3.3) 

tQ to 

yo 

{d^ + m^ + Lo{x)) p{x,y) = J dzI-P ix,z)piz,y) ■ 


(3.4) 


■^0 


These equafions are exacf, fhe approximafion is encoded info fhe used self energy, given by Eq. (p!^). 
Working ouf fhe self energy for fhe model of choice we can close Eqs. (3.3j^). Eor fhe one- 
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component scalar A<I>^/24 theory at three-loop level two diagrams contribute (at 0 = 0) [p 


r^(v,y) = -^F{x,y) (F\x,y)-^-p\x,y) 


X 


s {fj) = -—p{x,y) { F {x,y) - —p {x,y) 


1 


SoW = ^F{x,x) 


(3.5) 

(3.6) 

(3.7) 


The solution of the resulting closed set of equations is not accessible without numerical meth¬ 
ods. If the initial ensemble is spatially homogeneous, the two-point functions depend on five vari¬ 
ables {xo,yQ,x — f) only. 


The actual form of the equations of motion (3.3-3.4) depend on the type of the initial condi¬ 
tions. The most simple choice is to assume a known Gaussian density operator at initial time Iq, 


the presented form of Eqs. (3.3-3.4) reflect this choice. In this case the integrals start form Iq, and 


the numerical integration of Eqs. (3.3-3.4) is conceptually straightforward (but may be technically 
involved). 

The initial conditions for the second order differential equations can be given by setting F (vq = 
to,x',yo = to,y) and its first order time-derivatives. The p{xo = to,x\yo = to,y) function is zero 
identically, its derivative is constrained by Heisenberg‘s commutation relation: dx^pixo = to,x',yo = 
GJ) = 5(T-y). 

An initial particle distribution no{p) can be easily prescribed using Gaussian initial conditions: 


F{x^ ,y^-,p)\xP=y<^=,0 = 


«o(p) +1/2 


(On 


dxodyoF{x ,y \p)\x0^y0^tQ = [«o(p) + l/2]ft>p, 

= 0 , 


(3.8) 

(3.9) 
(3.10) 


with (Op = yjp^ +3^0 3Tq is the initial time physical mass. Here the propagators are Eourier 
transformed in space with respect to the relative coordinate x — y. 


4. The resummation scheme 

Eet us now review the resummations automatized by the 2PI formalism on the simple example 
of a scalar <J>^ theory. The considered approximation we define by giving the 2PI effective action 
explicitely in Eq. (2.6). The corresponding self energy equation (^) can be graphically expressed 


as seen in Eig. ^ The first two diagrams contribute at leading order to the self energy ^ The next 
two terms account for the next-to-leading order effects. 

By unrolling these equations and drawing the diagrams contributing to G in terms of the free 
(or initial time) propagator one reveals the elementary physical processes caught by this scheme. 
As shown in Eig. |3| a typical diagram is a ladder with finite number of rungs. Thinking of quasi¬ 
particles, each of these rungs represents an elementary scattering in the time dependent bath of 


*In the broken phase the external field (<I>) may carry inverse powers of the coupling, this will promote the <I>- 
dependent diagrams to lower orders. 
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Figure 2: Graphical representation of the 2PI equations of motion. The upper equation shows the self 
energy in the simple truncation of the scalar theory at three loop order. The thick lines stand for the 
self-consistent propagator (G), which is a subject of Dyson resummation (lower equation) using the free 
propagator (thin lines) and this self energy. 


- ) 0 « II ) • - 

^- 

time evolution 


Figure 3: Nonequilibrium interpretation of the ladder resummation. For any perturbative order there is 
a time when its contribution starts being relevant. A given order may account for some finite number of 
elementary collisions. After that number of elementary processes new orders will become of the same order. 
The 2PI ladder resummation considers all orders and solves the problem of secularity. 


other particles. One may associate a time scale T to these elementary processes. In a period of 
nx, n elementary collisions occur per particle on average. Beyond this period of time diagrams 
with more rungs start becoming relevant, hence higher perturbative orders are necessary. This 
phenomenon, called secularity, forbids the perturbative treatment of nonequilibrium fields, and 
calls for the resummation of this chain of ladders. This chain of ladders is what 2PI actually 
resums. 


If one now tries to calculate a four-point function out of £[(5] by cutting a line, one realizes 
soon that only one of the three possible channels are resummed, the Bose symmetry is broken. 
This is the reason we did not define fhe propagator as G, fhe solufion of fhe 2PI equations of 
mofion, buf used fhe IPI effecfive acfion insfead to introduce the proper two-point function W in 
Eq. dni)- The difference between G and ^ are diagrammatically exemplified in Fig. Q In many 
cases (see Eq. ( 2.12 )) fhey agree af 0 = 0, Iherefore we show a selecfion of ^-dependenf diagrams 
only. The channels missing from G can be included by solving a further (Bethe-Salpeter) equation 
]. This equation follows from the definition of W without additional theoretical input. 


26 


In Fig. 1^ we give its diagrammatic form, showing the NLO contribution only. One can show for 
the truncations obeying Eq. (2.12), that the Bethe-Salpeter equation is the equation of motion 
for at vanishing background, which is required to carry out the derivatives in 

Eq. (pd]). The four-point box contributes symmetrically in each channel to Eq. (|2.1[), this restores 
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2PI resummed IPI 


Figure 4: Graphs contributing to the 0-dependent part of the 2PI-resummed propagator at next-to-leading 
order. The 2PI equation of motion resums a ladder in one channel (left hand side). To the proper two- 
point function however, additional diagrams contribute (right hand side), that play an important role in 
the broken phase; they make sure that the Goldstone theorem is fulfilled. The missing diagrams can be 
calculated by solving the Bethe-Salpeter equation. 



Figure 5: Graphical representation of the Bethe-Salpeter equation. The thick lines stand for the solution 
of the propagator equation at vanishing background. 


the Bose symmetry [31, 22]. 

The topology of the resummed diagrams become more complicated for the 0{N) symmetric 
scalar theory in the large N expansion. As detailed in Refs. []_5, [I^, an additional resummation 


is required on the level of the 2P1 effective action. Then rint[0,G] is already an infinite series 
of diagrams Fig. ^ (left). By introducing a self-consistent vertex function (see Fig. ^ (right-top)) 
one may reduce the 2P1 effective action to a sum of finite number of 2P1 diagrams Fig. |6| (right- 
bottom). The concept of self-consistent vertex functions can be generalized and made systematic 
in the context of the n-particle irreducible effective actions [0^. 

This large-A resummation is always necessary if the occupation numbers are large (propor¬ 
tional to the inverse coupling). Any additional loop in Fig. ^ (left) carries additional powers of 
occupation numbers that compensate the coupling. This typically happens in preheating dynamics 
of the early universe 0- 

2P1 can be also used to derive Boltzmann equations. From the 2P1 effective action in Eq. ( p^ 
one arrives to a kinetic theory that accounts for the 2 —> 2 scattering of quasi-particles [|3^. The 
connection of the 2P1 equations with kinetic theories justifies the interpretation of the ladder re¬ 
summation above, and also makes clear what additional assumptions are made in such theories. 
Since 2P1 does not make these assumptions, it can be used to establish the range of validity of the 
following steps: 


1. The lower bound of the time integral in Eqs. (3.3-3.4) is extended to —oo. The equations are 
often formulated in Wigner transformed form. In order to make use of some mathematical 
identities to rewrite the equations in a handy form one usually extends the finite integrals in 
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Figure 6: Large-A^ resummation in the 0(N) scalar theory. The 2PI effective action to next-to-leading 
order in 1 /N consists of an infinite series of diagrams (left panel). These can be resummed by introducing a 
self-consistent coupling (big dot on the right panel) [|^ . 


the relative {x — y) coordinate to the range —oo.. . 00 . If one does this, the equations do not 
refer to an initial value problem any more. This introduces an approximation, that is justified 
from that time t on, when F{t,to) and dtp{tdo) become small compared to the equal time 
correlators. This usually occurs around damping time (tdamp), given by the imaginary part of 


the self energy [36]. 


2. Transport equations typically use gradient expansion to make the equations of motion local 
in time. It is controlled by the smoothness of the evolution, and becomes exact near equilib¬ 
rium. As expected, in weakly coupled theories subsequent orders of this expansion converge 
well. The intermediate time description of evolution can be substantially improved by going 


beyond the usually considered lowest order [136|]. 


3. The transport equation for p can be replaced by an assumed relationship between the F and 
p propagators. This nonequilibrium KMS condition parametrized by the particle distribution 
can be used as an ansatz in the F equation. It has been shown using 2PI equations that this 
relationship is dynamically established in damping time scale, well before the equilibration 
of the parametrizing distribution [37]. 


4. In the final step one replaces the spectral function (p) by a delta function, and introduces the 
quasi-particles. A detailed comparison of Boltzmann equation restricted to binary collisions 
with the 2PI dynamics based on Eq. (|2.6D can be found in Ref. [^. 


In summary, these steps are well justified in sufficiently weakly coupled theories, only after the 
damping time scale. This sets a clear limitation to transport and kinetic theories when they are 
applied to initial value problems [^]. 


5. Testing 2PI 

5.1 Convergence in equilibrium 

Since nonequilibrium description of fields requires non-perturbative treatment we were forced 
to use a resummation formalism. The same techniques can be applied in thermal equilibrium, 
too, where efficient formulations in Euclidean space-time become available. In contrast to the far- 
from-equilibrium case, there are various powerful approximation schemes known in thermal field 
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theory. A prominent approach in equilibrium high-temperature field theory is the so-called “hard- 
thermal-loop” resummation [pd]]. However, explicit calculations of thermodynamic quantities such 
as pressure or entropy typically reveal a poor convergence except for extremely small couplings. 
An important example for this behavior concerns high-temperature gauge theories. Recent strong 
efforts to improve the convergence aim at connecting to available lattice QCD results, for which 
high temperatures are difficult to achieve. In order to find improved approximation schemes it 
is important to note that the problem is not specific to gauge field theories. Indeed it has been 
documented in the literature in great detail that problems of convergence of perturbative approaches 
at high temperature can already be studied in simple scalar theories. For recent reviews in this 
context see Ref. |42]. 

A promising candidate for an improved convergence behavior is the loop or coupling expan¬ 
sion of the 2PI effective action. So far, thermodynamic quantities such as pressure or entropy 
have been mainly calculated to two-loop order. However, aspects of convergence can be sensefully 
discussed only beyond two-loop order since the one-loop high-temperature result corresponds to 
the free gas approximation. Efforts to calculate pressure to nontrivial order include so-called ap¬ 
proximately self-consistent approximations |^], as well as estimates based on further perturbative 
expansions in the coupling and a variational mass parameter [M]. These studies indicate already 
improved convergence properties. However, perturbatively motivated estimates as in Ref. suf¬ 
fer from the presence of nonrenormalizable, ultraviolet divergent contributions and the apparent 
breakdown of the approach beyond some value for the coupling. If one does not want to rely 
on these further assumptions, going beyond two-loop order requires the use of efficient numerical 
techniques. Such rigorous studies are important to get a decisive answer about the properties of 2PI 
expansions. As it turns out these problems appear as an artefact of the additional approximations 
employed and cannot be attributed to the 2PI loop expansion. 

In Ref. [ p^ ] we calculate the renormalized thermodynamics of a scalar <I>^ in the 2PI loop 
expansion to next-to-leading order. We solve the Schwinger-Dyson equation for the propagator, as 
well as the Bethe-Salpeter equation for the four-point function. By adding counterterms to the 2PI 
effective action we require the finiteness of the propagator following Refs. |^, ^]. In addition 
to this, we also renormalize the value of the effective action itself by an additional counterterm 
~ <I>^. As it has been shown in Ref. [22], all u-point functions are finite, and all counterterms are 
temperature independent. 

The bad convergence behavior of the perturbation theory is not seen in the 2PI results (Fig. 

To show this, we calculate the pressure to leading and next-to-leading order in the loop expansion, 
both in the 2PI and in the standard perturbative expansion. We could explore a broad range in the 
coupling, going very close to the triviality limit. 


5.2 Testing out of equilibrium 

Although the results from the 2PI formalism seem plausible, it has to be checked against other 
reliable nonequilibrium methods. If the exact evolution of the quantum theory were available, this 
test would be straightforward. 

In fact, the classical field theory approximation, provides the only way to solve the full non- 
perturbative dynamics of a field theory for arbitrary late times. In the classical limit the operator 
equations reduce to wave equations, which are convenient for numerical treatment. The classical 
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9r Qr 


Figure 7: Convergence of the 2PI resummation (left) and the perturbative series (right) in equilibrium. 
The scalar <I>^ theory has been renormalized at a temperature equal to the thermal mass Tq = mR{To), and 
also solved at T for which the pressure is shown. The subsequent orders of 2PI show a good convergence 
behavior even at higher couplings . 


approach became very popular both for heavy ion physics [45] and in cosmology [46], where 
the high occupation numbers keep the applications within the classical limit. The appearance of 
Rayleigh-Jeans divergences and the lack of genuine quantum effects, however, limit their use. 

Still, the classical approximation proved to be extremely useful for benchmarking the 2PI 
resummation in a nonequilibrium situation. The 2PI formalism is a generic field theoretical tool, 
that can be formulated both in quantum and in classical context. Aarts and Berges compared the 
classical and quantum dynamics of an 0(A) symmetric scalar field based on the large N expansion 
of the 2PI effective action, and solved the exact classical dynamics as well lo¬ 
in their comparison (see Fig. ^ for sufficiently high occupation numbers and large number of 
fields all the three approximations gave the same result. Even at lower occupancies and smaller N 
they found a good agreement between the classical 2PI and the exact classical results. This suggests 
that the selective resummation of the 2PI formalism includes the essential diagrams for irreversible 
dynamics. 

Future tests of 2PI dynamics might become available by a recent observation by Berges and 
Stamatescu [^]. They carried out a lattice simulation in Minkowski space-time by using a refor¬ 
mulation of stochastic quantization for the path integral [|4^. Although not much is known about 
the general convergence properties, this simulation technique is a promising non-perturbative ap¬ 
proach to out-of-equilibrium fields. A positive result of comparing 2PI with the lattice simulation 
would mark a breakthrough in nonequilibrium field theory [^]. 


6. Time scales of nonequilibrium evolution 

The 2PI equations of motion have been in use for several years, already. Most applications 

5 

mation but the energy density. All these applications confirmed that different quantities effectively 
thermalize on different time scales. Hence, a partial thermalization may be sufficient to support the 


21,M, 29, 30, 3^ i.e. the gradual loss of all initial infor- 


focused on thermalization [13,^, 15 
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Figure 8: Nonequilibrium evolution of the equal time propagator in a scalar 0(N) model for various 
momenta. Modes around pts were initially excited and their dynamics has been followed by solving the 
exact classical equations and the next-to-leading order 2PI equations both in the quantum theory and in the 
classical approximation. The agreement of the curves for = 10 shows that the 2PI techniques provide a 
very good approximation scheme even for late times. 0 - 

assumptions of thermal approaches in cosmology or in heavy ion physics. This has been pointed 
out in Ref. [||], where it was shown for a chiral quark-meson model that the prethermalization of 
important observables occurs on time scales dramatically shorter than the thermal equilibration 
time. As an example, Fig. ^ shows the nonequilibrium time evolution of fermion occupation num¬ 
ber for three different momentum modes in this model. The evolution is given for two different 
initial particle number distributions A and B shown in the insets, with same energy density. The 
vertical line marks the characteristic time scale ~ tdamp, after which the details about the initial 
distributions A or B are effectively lost. The following long-time behavior to thermal equilibrium 
is shown on a logarithmic scale in units of the scalar thermal mass m. 

In contrast to the very long time teq for complete thermal equilibration, prethermalization of 
the (average) equation of state sets in extremely rapidly on a time scale 

Ipt ^damp ^eq • (6.1) 

In Fig. we show the ratio of average pressure (trace over space-like components of the energy- 
momentum tensor) over energy density, w = p/e, as a function of time. One observes that an 
almost time-independent equation of state builds up very early, even though the system is still 
far from equilibrium! Here the prethermalization time ipt is of the order of the characteristic in¬ 
verse mass scale . This is a typical consequence of the loss of phase information by summing 
over oscillating functions with a sufficiently dense frequency spectrum. If the “temperature” (T), 
i.e. average kinetic energy per mode, sets the relevant scale one finds T tpt ~ 2 — 2.5 [Q]. For 
T > 400 — 500 MeV one obtains a very short prethermalization time tpt of somewhat less than 1 fm. 

This is consistent with very early hydrodynamic behavior, however, it is not sufficient as noted 
in Refs. Beyond the average equation of state, a crucial ingredient for the applicability 
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Figure 9: Fermion occupation number for three different momentum modes as a function of time in the 
chiral quark meson model of Ref. [plj]. 



Time [m'^] 

Figure 10: The ratio of average pressure over energy density w as a function of time. The inset shows 
the early stages for two different couplings h and demonstrates that the prethermalization time is rather 
independent of the interaction details. 


of hydrodynamics for collision experiments |§] is the approximate isotropy of the local pressure. 
More precisely, the diagonal (space-like) components of the local energy-momentum tensor have 
to be approximately equal. Of particular importance is the possible isotropization far from equi¬ 
librium. The relevant time scale for the early validity of hydrodynamics could then be set by 
the isotropization time. The analysis of scalar models lead to an isotropization time given by the 
comparably long characteristic damping time ~ tdamp P^]- In gauge theories, however, there is 
a weak-coupling mechanism for faster isotropization identified as plasma instabilities 0i01- 
Whether this can explain the observations or whether they suggest that we have to deal with some 
new form of a “strongly coupled Quark Gluon Plasma” is an important open question. 
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